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The spectra and electroweak decay properties of light 
mesons are analyzed within the framework of the instanta- 
neous Bethe-Salpeter equation. The interaction kernel com- 
prises alternative spin-structures for a parameterization of 
confinement and a residual quark-antiquark interaction based 
on instanton effects. It is shown that only with a vector con- 
finement the parameters can be chosen such as to yield an ex- 
cellent description of the light pseudoscalar and vector mesons 
including weak and two photon decays. However it is found 
that it is not possible to reconcile this with the Regge be- 
havior of higher lying meson states with the same parameter 
set. 



I. INTRODUCTION 

In a previous paper [Q it has been shown that the 
Bcthc-Salpetcr(BS) equation with an instantaneous in- 
teraction (Salpeter equation) provides a suitable frame- 
work for relativistic quark models of light mesons. The 
use of the Salpeter equation has several advantages com- 
pared to other treatments, i.e. 

• The relativistic kinematics of the quarks is treated 
correctly; 

• The amplitudes have the correct relativistic nor- 
malization; 

• The lower component $ of the Salpeter ampli- 
tude is determined dynamically. 

On the other hand the practical advantages of a nonrel- 
ativistic treatment are also present, i.e. 

• The Salpeter equation can be formulated as an 
eigenvalue problem Hip = Mip for the mass M of 
the bound state. 

• One can define a (not positive definite) scalar prod- 
uct (ipal^b) for t ne Salpeter amplitudes. 

• The Salpeter operator Tt is selfadjoint with respect 
to this scalar product. 

In [jlj we have developed a numerical scheme which en- 
ables the calculation of meson mass spectra and the cor- 
responding BS-amplitudes. Furthermore we have given 
formulas to compute some important electroweak meson 
decay widths. 



In the present paper we shall apply this method to an 
explicit quark model for light mesons. Our main concern 
in this context will be whether a realistic description of 
deeply bound states like the pion is compatible with a 
reasonable description of confinement. 

Let us first give a list of the main features of light 
mesons that will be considered in the following: 

• the low masses of tt and K 

• the weak decay constants / w and fx 

• the decays 7r°, rj, v( — > 27 

• the masses and the flavor mixing coefficients of 77 
and rj 

• the masses and the leptonic decay widths for the 
p, to and (f> mesons 



• the Regge behavior M 2 ~ J 

As far as we know there presently is no model that can 
describe all these features in a consistent way. On the one 
hand there is the nonrelativistic quark model that gives 
a reasonable description of the mass spectra , but that 
completely fails in describing the decay width s of the 
deeply bound states like the pion (see Sec. IV D). On the 
other hand there are models like the Nambu Jona-Lasinio 
model ||-[5| that are based on the chiral symmetry of 
QCD for vanishing current quark masses. This model 
leads to a good description for the tt, K, rj and rj mesons, 
but higher angular momenta states or radial excitations 
cannot be described since confinement is ignored. 

An attempt to arrive at a more complete description 
based on the Salpeter equation has recently been given 
by J.F.Lagae f|. His results show the difficulty of finding 
a suitable ansatz for the confining interaction kernel: he 
rejects the hypothesis of an instantaneous scalar confin- 
ing kernel. 

The model we present in the following is based on a 
linear scalar or alternatively vector confining kernel com- 
bined with an effective interaction computed by 't Hooft 
from instanton effects in QCD |?]-[|]. A nonrelativistic 
version of this interaction has already lead to good re- 
sults for the meson and baryon mass spectra We 
therefore feel encouraged to test this ansatz in the rela- 
tivistic Bethe-Salpeter framework. 

The paper is organized as follows: In Sec.[n]we briefly 
summarize the main features of the Salpeter equation 
as given in ref. (!]]. The explicit form of the confining 
BS-kernel and of the 't Hooft kernel is given in Sec. III. 
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In Sec. IV we present the calculated meson mass spectra 



and obtain the pion and kaon decay constants / ff , fx, the 
decay width into two photons for the 7r°, rj and T)' mesons 
and the leptonic widths for the vector mesons p, lu and 
</>. We will compare these d ecay w idths to corresponding 
nonrclativistic results in Sec. IV D using the wave function 
of ref. which shows the impressive improvement due 
to the relativistic treatment of the quarks compared to 
the nonrelativistic potential model. Finally we give some 
concluding remarks in Sec.|v|. 



II. THE SALPETER EQUATION 

For an instantaneous BS-kernel and free quark propa- 
gators the p° integrals in the BS-equation can be com- 
puted analytically in the rest frame of the bound state 
with mass M. The result is known as the Salpeter equa- 
tion which reads 



d 3 p> Ar(p) 7° PW) *(P')] 7° A+(-p) 



(27r) 3 M + W\ + LU 2 

d 3 p> A+ (p) 7 ° l(V(p, p') Hp')} 7 ° A 2 ~ (-P) ( 

(2tt) 3 M -UJX-0J2 



with uji = ^Jp 2 + mf, Af (jp) = (u, ± Hi(p))/(2u)i) and 

Hiip) = 7°(7P + m i)- 

This equation can be reformulated into an eigenvalue 
problem for the bound state mass M. With the defini- 
tions 

V>($ :=<!>($ 7° (2) 
[W(p,p') W)] := 7° [V(p,p>) (3) 

the Salpeter equation can now be written as 

Hip = Mip with 
Hip = H x ip - iPH 2 

{A+[^')]A 2 --A7[^')]A 2 + } (4) 

with all quantities depending on p or p' (indicated with 
a prime) and the notation J = J d 3 p' /(2n) 3 . 

One can now define a (not positive definite) scalar 
product for amplitudes ipi = $17° and 1P2 = $27° as 

(iP, \ip 2 ) = J tr (ip\ A+ ip 2 A 2 - - 4 Ar V2 A+) (5) 

The normalization condition for solutions of the Salpeter 
equation is then given as (ip\ip) — {2tt) 2 2M. If one con- 
siders amplitudes satisfying A^ipA^ = A^ipA^ = one 
easily finds 



(1P1 \ Hip 2 ) = I (w!+co 2 )tr (ip\ip 2 ^ 

■x (i>\ WiP' 2 ) (6) 



For most kernels of physical interests H is self adjoint with 
respect to the scalar product given above which implies 
that bound state masses M are real for eigen functions 
ip with nonzero norm and that amplitudes ipi and ip 2 
corresponding to different eigenvalues M% ^ M 2 are or- 
thogonal, i.e. (ipi I ip 2 ) = 0. 



III. THE BETHE-SALPETER KERNEL 
A. Confinement 

Up to now the confining interaction of QCD is only 
known in the static limit of heavy quarks. In this limit it 
has been shown that the static potential between 

quarks is of the form Vc (r) = a c + b c r + W where W de- 
notes the relativistic corrections of the order p 2 /m 2 . As 
stated by Gromes (l(]] it is still an open question whether 
an additional \ jr term should also be included into the 
static confining potential. Usually one concludes from 
the sign of the spin-orbit coupling term in W that the 
confining qg-interaction behaves like a Lorentz scalar. 

The problem for the case of light quarks is that up to 
now there is no unambiguous extension of the confining 
potential beyond the static limit. Especially there is no 
prescription on how to extend it to a noninstantaneous 
form. Naive noninstantaneous extensions fail as has been 
shown by S.N.Biswas et al. (T^| for the harmonic oscilla- 
tor BS-kernel V(x) = ~bx 2 = b(x 2 - (a; ) 2 ) that yields 
only a continuous spectrum. Similar results are to be 
expected for other kernels like 1/q 4 - Because of these 
difficulties the only way we see at the moment is to pa- 
rameterize confinement as an instantaneous interaction 
kernel. 

The sign of the LS-term in the static limit would be 
compatible with a scalar confinement kernel. However 
some authors [p 13 have shown that the linear Regge be- 
havior M 2 ~ J is lost for this choice, since bound state 
masses come out too small for higher angular momenta. 
This is due to the relativistic corrections to the static 
potential and becomes more problematic with decreasing 
quark masses. For a vector confining kernel this problem 
almost disappears, but in the static limit a vector kernel 
leads to an LS-term which has the wrong sign. To our 
knowledge there presently is no convincing parameteriza- 
tion for the confining kernel that exhibits both features, 
i.e. leads to linear Regge trajectories and yields a spin- 
orbit term with the correct sign. 

In the following we will analyze both spin structures 
for the confinement, i.e. a scalar 1 (g> 1 and a vector 
7° ® 7o interaction. In the rest frame of the bound state 
the corresponding BS-kernels in eq.(|l|) are parameterized 
as 

[V^(p,p r ) Hp')] = -VcW- P'f) 7° W) 7° (7) 

[vS(p,p')Hp')] = v c {{p-p'?)W) (8) 
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where Vc is a scalar function which has the fourier trans- 
form V(j(r) — a c + b c r. 



B. 't Hooft interaction 

1. The 't Hooft lagrangian 

It is already known from nonrelativistic potential mod- 
els that the masses of the scalar and pseudoscalar mesons 
7r, K, 77, 7/ cannot be described with a confining poten- 
tial alone. The usual extension would be to add another 
contribution to the interaction kernel that comes from 
a One Gluon Exchange (OGE). This works quite well 
for heavy quarkonia Jl3|,|l4| , but it is highly questionable 
for light mesons where perturbation theory cannot be ex- 
pected to work. Since OGE leads to a flavor independent 
interaction kernel one thus obtains degenerate it and n 
mesons in clear contradiction to the experimental mass 
values m?; = 140 MeV and m v = 549 MeV. In order to 
cure this discrepancy one would have to take into account 
higher order diagrams. 

However there is another QCD based candidate for a 
residual qq-interaction computed by 't Hooft and others 
from instanton effects (7]-|j| which leads to good results 
for meson and baryon mass spectra within a nonrelativis- 
tic potential model 

Instantons are special solutions of the classical non- 
abelian Yang-Mills equations in Euclidian space. They 
are peaked both in space and imaginary time having a 
finite extension p. Since they cannot be deformed con- 
tinuously into classical solutions corresponding to gluon 
fields they lead to an effective interaction between quarks 
that is not covered by perturbative gluon diagrams. This 
interaction leads to spontaneous breaking of chiral sym- 
metry as can be seen by normal ordering of the under- 
lying Lagrangian (see appendix [A|) . The normal ordered 
Lagrangian takes the form 

3 

C = k + ■ (%7 M <^& - nij qq) : +A£(2) + A£(3) 

(9) 

where k is an inessential constant that renormalizes the 
vacuum energy. A£(2) and A£(3) are two and three 



body terms and nij 



+ Am, is the effective con- 



stituent quark mass. In the following we will consider 
these terms. 



2. The constituent quark mass 



The contribution Am, to the effective constituent 



quark mass m. 



Am, is given by 



d (p) 4 2 „ / 2 2 



Am„ = J dp—^---K p \m n p- -n p (q n q„) 




d (p) 4 2 3 / 2 2 .; 



Am s = j dp^-^p A [mlp--n'p- i (q n q n )j (10) 


where the instanton density for three colors and three 
flavors reads |ll| 



do(p) = (3.63 ■ 10~ 3 ) 
with 



Sir 2 



Sir 2 



9>(p)j exp \rm: 1111 



87rM =91nf-^ + ^minf-^) (12) 



9 2 {p) 



A QCnP/ 9 K^QcnP 



within two loop accuracy [Q. Here A QCD is the QCD 
scale parameter. The integration over the instanton size 
p has to be carried out up to a cutoff value p c where the 
In In term coming from the two loop correction is still 
small compared to the In term. 

Equations ( [To| ) are usually called Gap equations. They 
describe the generation of a dynamical quark mass due 
to the interaction with the negative Dirac see. In our 
model the constituent quark masses nij will be used as 
free parameters that are fitted to the experimental data. 



It will be checked in the end (see Sec. IV E) whether the 



obtained quark masses are compatible with common val- 
ues for the quark condensates present in Ato^ assuming 
that the confinement interaction does not contribute es- 
sentially to the process of chiral symmetry breaking. 



3. The two body interaction 
The two body term reads (see appendix |A| ) 



A£(2) =-^EEE #eff(*) £ ikl£zmn 



i kl mn 



{ : qlqj [70-70 + 7o75-7o7s] (27f + P 6 C ) q m q n : } (13) 
where the effective coupling constants are given as 

g eS (i) = Jdp^l (%y) 2 (m°p- |*V<ft«: 



(14) 



and the tensor notation 



ij k.l 
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has been used for Dirac and color indices. This repre- 
sentation explicitly shows the antisymmetric flavor de- 
pendence of the interaction. The color sextett and an- 
titriplett projection matrices are given by 



V§ = \(l C 



u c ) = 2 -i c + 1 -x-x 



so that 



2Vk 



T° = -(31 c -IT C ) = -1 C --X-X 
6 2 v j 3 4 



(16) 
(17) 

(18) 



where H c is a color exchange matrix defined as H c 
da Sjk and A a (a = 1, . . 



are the color matrices. Just 



like the quark masses also the coupling constants will be 
treated as free parameters in our model. 



4- The three body interaction 

After a lengthy calculation the three body force can 
finally be written in the form 



A£(3) 



27 



.0). 



80 5 eff1 
[70 • 7o • 7o + 7o75 • 7o 75 • 7o 
+7o75 • 7o • 7o75 + 7o • 7o7s ■ 7o7s] 
V[(2V C + 5Vg)qqq :} 



(19) 



where Vf is the projector onto a three-particle flavor 
singulett state, Vxq and V§ are projectors onto the color 
decuplett and the color octett. The effective three-body 
coupling constant is given by 



7 (3) 

? cff 



Pc , MP) 

P 



b 2 3 
-7T p 



(20) 



Obviously this three body force does not contribute to 
qq-states and to colorfree qqq-st&tcs. 



5. The t Hooft kernel 

In order to distinguish different indices in the following 
we will use the notation Sj for Dirac indices, Cj for color 
indices and /, for flavor indices. The vertex correspond- 
ing to AC2 is shown in Fig.[l| The vertex reads 



(H) G /3/4,/l/2 (l S3 Slls 4 S 2 +Js 3Sl ls 4 S 2 ) 

(4/3 lc 3 ci lc 4 c 2 1/4 X C3Cl X C4 c 2 ) 
with the definition 

G /3/4,/l/ 2 - -9eff(/5) e /5/ 3 /4 C /5/l/2 

8 h 



(21) 



(22) 



For the qq — > qq amplitude we have to consider the two 
diagrams given in Fig.||. In a meson the quark and the 
antiquark are in a color singulett state represented by the 
matrix x? lC2 = ^cic 2 /V3- The color matrix elements for 
the two vertices are then given by 

1 

2 



C1C2 C3C4 



c 

C3C4 



) 9 (3 lc 3 ci lc 2 c 4 nc 3 c 2 ,cic 4 ) X ClC2 — (23) 



C 

C3C4 



1 



n 



c 3 c 2 ,c 4 c 



1 )Xc lC2 = 4 ( 24 ) 



so that only the second vertex in FigJ^ contributes to 
the interaction kernel. The effective 't Hooft interaction 
vertex between qq color singulett states is then given by 



(-4i)G /2/3 ,/ l/4 (l f 



,1. 



si 7s 3 s 4 7s 2 si ) 



(25) 



Since we assume 5'?7(2)-flavor invariance of the interac- 
tion we set 



9 



fs e ff( s ) 



9 



(26) 



where s stands for strange and n stands for nonstrange 
(u,d) flavor. The results for Gf 2 f 3 j 1 f i are given in Tab.|. 
Tab-H shows the matrix elements of Gf 2 f 3 j 1 f i for the 
flavor functions 



7T° = (uu - dd) /a/2 
r/ n = (uu + dd) I a/2 



'Is 



(27) 
(28) 
(29) 



From the vertex we can extract the lowest order contri- 
bution of the 't Hooft interaction to the BS-kernel as 

[V T (p,p>) $(p')] hh = 4 J2 GhfiMi 
flfi 

[1 tr ($ /{/ , (p')) + 7 5 tr (p 1 ) 7 5 ) ] (30) 

As shown in the appendix this interaction only acts on 
the scalar and pseudoscalar mesons 

j-k p = ±^ For the 

pseudoscalar mesons it is attractive for the it with a cou- 
pling constant g and for the K with a coupling constant 
g' . For the r\ and 77' the interaction leads to mixing of 
nonstrange and strange flavor amplitudes. The effective 
sign of the interaction is reversed for the scalar mesons 
thus being repulsive for the ciq. Note that in the nonrel- 
ativistic limit the 't Hooft interaction only acts on the 
pseudoscalar states J 7rp = 0~. The 't Hooft kernel as it 
stands represents a pointlike interaction that has to be 
regularized. Following ref. j^] we do this by multiplying 
the kernel with a regularizing Gaussian function 



Vreg(<?) 



(31) 



with q = p — p' and q = \q\. In coordinate space this 
choice corresponds to replacing the 8(f) function by 



which introduces a finite effective range A. 



(32) 
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IV. RESULTS AND DISCUSSION 

A. Models and Parameters 

The main concern of our work was to see whether we 
can obtain a consistent description of a) the masses and 
decays of the low lying pseudoscalar and vector mesons 
and b) confinement reflected e.g. by the Regge trajecto- 
ries. 

For this purpose we investigate two different models 
of the confinement kernel: 1) a vector 7 ® 7 - and 2) a 
scalar 1 ® 1-structure. 

The parameters used are the nonstrange and strange 
quark masses m„ and m s , the offset a c and slope b c of 
the confinement interaction, the two coupling constants 
g, g' and the effective range A of the residual instanton 
induced interaction. So the total number of parameters 
amounts to seven. 

We used two sets of parameters in the vector confine- 
ment case: Model VI was tuned to reproduce the masses 
and decays of the low lying mesons. We therefore used 
a small nonstrange quark mass m n , as the correct de- 
scription of the pseudoscalar decays depends essentially 
on this quantity. Given this mass we had to take a mod- 
erate confinement slope to reproduce the decays of the 
vectormesons. The offset a c was fixed by the p-mass and 
m s by the if*-mass. Finally g and g' were fixed by the 
masses of the pseudoscalars tt, 77 and K. In Model V2 we 
used a larger nonstrange mass of m n of about 1 /3 of the 
nucleon mass, which is a value common to Nonrelativis- 
tic Quark Models. The aim of this parameter set was to 
obtain a good description of the Regge trajectories and 
the higher lying resonances. 

Finally in Model S we investigated a scalar confine- 
ment provided with the same quark mass m n as in V2. 
For this kernel it turned out that our method of solv- 
ing this type of BS equation works reasonably for higher 
angular momenta only if the fraction of the confinement 
slope b c and the quark mass m q is sufficiently small. The 
parameters for the three model are listed in Tab. HI. 



state with imaginary mass and zero norm. The physical 
interpretation of this phenomenon however is not clear. 
For a scalar confinement this effect is compensated by 
the negative spin orbit splitting. In the following we will 
discuss the differences of the three parameter sets. 

Since in model VI we used a small confinement slope to 
reproduce the decays of the vectormesons, the calculated 
Regge trajectory is too flat (see Fig.|[). The spin orbit 
splitting between the 1 ++ and 2 ++ mesons in our model 
is purely due to the confining interaction. In both models 
with vector kernel it is of order of 200 MeV exceeding 
the experimental mass difference, which is in fact rather 
small (see TabQ. 

In model V2 with large quark masses and large confine- 
ment slope b c we obtain a good description of the masses 
of all mesons comparable to the results in nonrelativistic 
calculations. However b c has to be much larger in the 
BS framework, as obviously the kinetic energy is overes- 
timated in nonrelativistic calculations. The Regge tra- 
jectories representing the confinement property are well 
reproduced (FigJ|). For the spin orbit splitting the same 
remarks as for VI apply. 

Finally model S shows a reasonable description of the 
ground states of the pseudoscalar and vector mesons. For 
the states with large angular momentum the method of 
solving the BS equation by a basis expansion does not 
lead to convergent solutions with positive norm. With 
increasing dimension of the basis the smallest positive 
eigenvalue decreases until it becomes imaginary (a simi- 
lar problem arises if one studies a 7^ (8 7^ interaction) . 
We therefore cannot explicitly exclude a scalar confining 
potential, mainly because with the present parameters 
we were not able to find stationary solutions for higher 
angular momenta. For these reasons we also omitted in 
the Regge plot the state with angular momentum j=4. 
The /1 meson mass is larger than the mass of the /2 me- 
son. Taken together with the results for the vector case 
this indicates that the spin orbit splitting can only be 
explained by a mixture of scalar and vector type interac- 
tion. 



B. Mass spectra 



C. Decay Observables 



The quality of the mass spectra is different for the three 
models, as different priorities led us to the parameters. 
Common to all three models is an overall agreement of 
the masses of the pseudoscalars 7r, X, 77,7/ and the vec- 
tor mesons p,lu,4>, K* . The spectra for these meson are 
compared to experiment in Figs. The 't Hooft inter- 
action leads to the correct splitting of 7r, 77 and 77' mass. In 
contrast with experiment, where the a and / are nearly 
degenerate, we obtain a large splitting due to the instan- 
ton induced interaction of several hundred MeV (compare 
Tab. IV). For a vector type kernel with positive spin or- 



bit splitting this leads an enormous attraction: although 
there is a scalar state at roughly lGeV, we also find a 



In this section we will discuss the influence of the pa- 
rameters on the decay observables of the pseudoscalar 
and vector mesons. 

The parameters of Model VI have been chosen in order 
to give a good description of the masses and the decays 
of the pseudoscalar and vector ground states. As can be 
seen in Tab.[v], we obtain an almost quantitative agree- 
ment. As an important result we consider the fact that 
the 7r and 77, which usually are interpreted as Goldstone 
particles of the of chiral symmetry breaking, also may be 
understood as bound states of quark and antiquark, al- 
beit with relatively small quark masses. This is reflected 
in the simultaneous agreement we obtain for the pion de- 
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cay constant /„. and the decay width tt° — > 77, which 
in the alternative Goldstone picture are related to the 
Adler Bell Jackiv anomaly ]l7j] . It becomes clear that 
a relativistic treatment of the decay formula and of the 
normalization are important for a correct understanding 
of the pion as deeply bound quark antiquark state. This 
can also be seen in Fig.|] showing the upper and lower 
component $ ++ , $ of the pion amplitude. In contrast 
to the wave function for the p meson (Fig]|) the upper 
pion amplitude is only about 10% larger than the lower 
one in the region of small relative momenta. For higher 
momenta the amplitudes even become equal. This fact 
obviously leads to important cancellations for the nor- 
malization and for the decay constant. This was already 
emphasized in an early quark model [^s| , where the ef- 
fects of different estimates for the relation between the 
upper and lower component on weak decay constants for 
pion and kaon were analyzed. For a comp ariso n with non- 
relativistic decay formulas compare Sec. IV D . Although 
the pion amplitude has significant contributions up to 
momenta of about 4GeV/c, the main part e.g. of the 
integral for the ir° decay width comes from momenta of 
about 150 MeV/c (as this is the scale of nonstrange and 
pion mass). 

For the 77 — > 77 decay we also find excellent agreement 
with experiment, whereas the process for the 77' is under- 
estimated. The results depend strongly on the correct 
nn- ss mixing, as e.g. for the 77 we obtain a negative 
interference. The mixing (due to the instanton induced 
interaction) can be compared to a simple model given by 
Rosner |l{J . The physical mesons are expanded in a basis 
of three states \N >= 1 / \/2 \uu + dd > , \S >= \ss > and 
\G >= \Gluonium >: 



\<r) >= X n \N > +Y V \S > 



Z rt \G> 



(33) 



The coefficients may be estimated from electromagnetic 
transitions p(J . We compared the results for the absolute 
values of X and Y in Tab. VI with the contributions of 



the nonstrange and strange part of the amplitude to the 
relativistic norm. The results agree well in the case of 
the 77, but not for the 77'. Experimental results indicate a 
larger gluonic component for the 77', which could modify 
the results. 

The leptonic decay widths for the vector mesons are 
also in good agreement with the data. This is essentially 
due to the small confinement slope, which determines the 
size of these mesons. We conclude that a consistent de- 
scription of all the ground state pseudoscalar and vector 
mesons is possible in this framework. However, more ob- 
servables like electromagnetic transitions or the pion form 
factor will be calculated in the future to substantiate this 
statement. 

The agreement with experimental data for model V2, 
which was designed to reproduce the higher resonances 
and Regge trajectories, is only of a qualitative character 
(Tab.0). Decay constants and photon decay widths dis- 
agree by about 50%, which is essentially due to the large 



quark mass. The leptonic decay widths of the vector 
mesons are overestimated due to the steep confinement, 
which enlarges the amplitudes at the origin in coordinate 
space. 

Comparing the two parameter sets for vector confine- 
ment we find: VI with light quark masses gives a quanti- 
tative description for the vector and pseudoscalar ground 
states, but only a qualitative picture of the Regge be- 
havior. For V2 with large quark masses the situation is 
opposite. This might be an indication that for large dis- 
tances the mass of the quarks effectively should increase 
due to some additional mass of a string. 

For model S we find good agreement for the vector 
mesons, but not as good for the pseudoscalars. Although 
it is possible to describe the latter with a smaller quark 
mass, we could not obtain a quantitative adjustment for 
both _ and 1~ with a scalar confining kernel. 



D. Comparison with nonrelativistic results 

In order to estimate the relevance of relativistic ef- 
fects in our model it is useful to compare the results with 
the corresponding ones computed in the nonrelativistic 
quark model. In the nonrelativistic limit the Salpeter 
equation reduces to the usual Schrodinger equation with 
the hamiltonian given by 



H = a c + b c 



\G 



b\r) 



(34) 



where the ^-function again has to be regularized accord- 
ing to eq. ([32"|) . The mass spectra for this hamiltonian 
have already been investigated in an earlier work being 
in good agreement with experiment. In the following we 
will use the wave functions i\>{f) of || to calculate the de- 
cay observables using the well known formulas ||l4| , f2l| , p2"f 



r(i- 
r(o~ 



/ W = ^J |^ ( 0)| 



'M 
16 it a 2 e 2 



i+r) = 

-> 77) 



M 2 
Uira 2 g 4 



(35) 
(36) 

(37) 
(38) 



Here M is the experimental meson mass, m q is the quark 
mass, a = 1/137 and e q gives the quark charge in units 
of the proton charge according to the flavor composition 
of the meson. The results for these decays are given in 
Tab.|v|. One finds that leptonic decays can already be 
described reasonably in a nonrelativistic framework. On 
the other hand the weak decay constants and especially 
the two photon decay widths are far away from the ex- 
perimental data. This discrepancy cannot be cured by 
changing the parameters within reasonable limits. We 
therefore conclude that a relativistic treatment is essen- 
tial for the pseudoscalar mesons. 
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E. Discussion of the gap equations 

Due to the process of chiral symmetry breaking the 
't Hooft interaction leads to relations for the effective 
constituent quark masses m n , m s and the coupling con- 
stants g, g' as given in eqs.(flO|). In order to check if 
these relations are qualitatively compatible with the fit- 
ted parameter sets we use A QCD = 200 MeV, m° — 
9MeV, to" = 150 MeV, («n?n) = (-225 MeV) 3 and 
(q s q s } = 0.8(g n g„) (compare and plot m n , m S) g 

and g' as functions of the instanton size cutoff p c as shown 
in Figs.|j| 

Because of the delicate dependence on the condensate 
values and due to the regularization procedure in the 't 
Hooft kernel one should not expect quantitative agree- 
ment with the fitted parameter sets. 

For p c = 0.408 fm one finds e.g. m n = 170 MeV, 
m s = 270 MeV, g = 79 MeV fm 3 and g' = 
58 MeV fm 3 . Despite of the strange quark mass which 
comes out too small the other parameters are quite close 
to the values of parameter set VI. Note that p c is almost 
equal to the effective range A of the 't Hooft interaction. 
Furthermore we find that g' is smaller than g for all val- 
ues of p c as is the case for all three parameter sets. 



V. SUMMARY AND CONCLUSION 

Within the framework of the instantaneous Bethc 
Salpeter equation we investigated different models for 
mesons as bound states of quark and antiquark. We con- 
sider it to be the main result of this work that the masses, 
weak decay constants and two photon widths of the light 
pseudoscalar mesons (ir, r\, K) can be understood quan- 
titatively in terms of a qq description, alternatively to 
the Goldstone picture. With the same parameters also 
the masses and leptonic decays of the vector mesons can 
be reproduced. 

The pseudoscalar mesons are dominantly affected by 
an instanton induced interaction, which apart from the 
7T, r\ splitting gives the correct nn— , ss— mixing for the 
77 meson. The self energy corrections and coupling con- 
stants due to the resulting chiral symmetry breaking are 
compatible with the parameters we used for the quark 
antiquark interaction. In contrast to nonrelativistic cal- 
culations instanton effects appear also in the scalar sec- 
tor and lead to a isoscalar state with imaginary mass and 
zero norm. 

Concerning the nature of the confinement kernel we 
find that a vector type interaction can reproduce the 
Regge trajectories, although with a larger quark mass 
than the one needed to describe the lowest lying mesons. 
This might be an indication that due to string effects the 
quark mass should increase with distance. 

We were only able to find reliable solutions with a 
scalar kernel for relatively large quark masses and weak 



confinement. With these parameters, however, one can- 
not reproduce quantitatively the ground state mesons or 
the Regge trajectories. Nevertheless the spin orbit split- 
ting indicates the existence of a scalar component in the 
interaction in order to cancel the large mass differences 
coming from the vector structure. 

Acknowledgments: We are thankful to M.Fuchs for 
constant help in numerical questions. 



APPENDIX A: THE 'T HOOFT INTERACTION 

1. The Lagrangian 

As shown by Shifman, Vainshtein and Zakharov || the 
contribution of an instanton-antiinstanton configuration 
to the effective quark Lagrangian for three quark flavors 
u, d, s is given by 



AC = dp 



32 V 3 
3 



ttV 



[{qiR\ a qiL){q2R\ a q2L} 
j(qiRa^X a q 1L )(q 2 Ra^X a q 2L )} 



mlp- g7r 2 p 3 (9 3it (j3 L ) 



ttV d abc {q 1R a^X a q 1L ) 



9_ 

+ 40 V3 

(q2R<J^X b q 2L )(q 3R X c q 3L ) 

+cycl. perm, of (123) 



"25g« I 3^ P ) / (giil^A qiL) 



(Q2RcrJX q2L){q3R^ T ^X c q 3L ) 



320 \3 

(q~3RX c q 3 L 



7T 2 p 3 ) d abc (q 1R X a q 1L )(q 2R X»q 2L ) 

^(R^L)\ (Al) 



with a^u := bv,7„] /2 and q lL := |(1 + 75)*, ffifl := 
5(1 — 75)9i being the projections of the quark Dirac op- 
erators qi onto left and right handed components. Fur- 
thermore i = 1,2,3 = u,d,s denotes the flavor degrees 
of freedom, to° the corresponding current quark masses, 
A a (a = 1, ... ,8) are the color matrices and f abc , d abc 
are the standard SU(3) structure constants defined by 
the commutator [A a , A b ] _ = 2i f abc X c and the anticom- 
mutator [A a ,A b ]+ = f<5 ab + 2d abc X c . The 't Hooft in- 
teraction leads to spontaneous chiral symmetry break- 
ing as can be seen by normal ordering C = Co + AC 



7 



with respect to the physical QCD vacuum, where Co = 
2j=i (^Qjl^^tJ.Qj ~ m j Ql) i s the free quark Lagrangian. 
Using the Wick theorem one finally obtains eq.(0). 



2. The interaction between two quarks 

After normal ordering the two body interaction term 
is given by 

A£(2) = 5 e ff( 3 ){ : (3iR9li92R?2i) : 

3 r 

+ 32 : |^ii?A a qiL)fei?.A a g2 L) 

3 1 
-^{qiRO l ivX l qiL){q2R(J^X 1 q2L)\ i ■ 

+{R <-> L)\ + cycl. perm, of (123) (A2) 



with 3 e ff(i) given in eq.(14). One can transform the two 
body force into a more transparent form using the no- 
tation Eijk,i — U,d,s with s u d s = 1- Insert qiL,R — 
(1 ± 75)/2 qi and use the relations 



-4 + 75 £*- 75 £ l ) 

(A3) 



with E = diag(a, a) where the notation (A ■ B)^ ^ — 
Aik Bji has been used. Furthermore use 



211 s - I s 



A a • A a = 2 n 



-T 



(A4) 



with IT , IL F and II being exchange operators in spin, 
flavor and color defined as Uy^; = 8a Sjk On the anti- 
symmetric tensors one has H H F H = — 1 which can be 
used to eliminate the spin dependence leading to eq.(13). 



APPENDIX B: NUMERICAL SOLUTION OF 
THE SALPETER EQUATION 

1. General framework 

In this section we briefly sketch the numerical method 
outlined in re f. [pj . With the standard Dirac representa- 
tion (see e.g. p4fl) $ is a 4x4-matrix in spinor space that 
can be written in block matrix form as 



$($++,$—) 



(Bl) 



where each component is a 2x2-matrix. The definition 
of the bilinear function <t is motivated by the fact that 
only two of these four amplitudes are independent since 
the Salpeter equation implies 



C 2 S$ 



(B2) 
Ci = 



where we use the shorthand notation s = ap, 
w l /(wim 2 +w 2 mi). 

The transformation properties of the Salpeter ampli- 
tude imply that we can decompose <I> ++ , as 



-(f) = 
-(P) = 



L S 

E 

LS 



n[- s \ P ) [Y L {n p )® Vs \ ] 



(B3) 



with the spin matrix (</?Sg icr 2 ) mm ' = (1/2 m l/2m'\S q). 
L, S have to respect the usual constraints coming from 
parity and charge conjugation. Let 

Et (p) = R n%u (P) [Y Li (fip) ® <ps t ] J Uj (B4) 

be a complete set of 2x2 basis functions orthonormal 
with respect to the scalar product given by (Ei \Ej) = 



used to expand $ 



= Sa. These basis functions can be 
with corresponding (usu- 



+ and $ 



ally real) expansion coefficients and a\ . Inserting 
these expansions into $ gives 



(B5) 



with e\ +) = $ (Ei,0) 7 ° and e\ ' = $ (0, EA 7 °. Defining 
the matrix elements (which are real in our case) Hf? = 



(-) 



\He 



(*')\ 



and N, 



the Salpeter equation 



TLip = Mip can now be written in the form of a matrix 
equation as 

'H++ H+~\ /a<+>\ _ „ (N++ 
M + - H++ L(-)i -JV+ 



(B6) 

We solve this equation within a finite basis i < i ma x = 10 
and use the variational principle SM = looking for sta- 
tionary points of M as a function of /?, where /3 is a vari- 
ational parameter with the physical dimension MeV~ x 
that sets the absolute scale of the basis functions as 
Ei(p) — /3 3 ^ 2 £"f _1 (p/3). For the basis functions in mo- 
mentum space the following functions have been used: 



,-v/ 2 



(B7) 



with N n L being a normalization coefficient, y = p/3 and 
L 2 / J+2 (y) being a Laguerre polynomial. The fouriertrans- 
formed basis states also can be given in analytical form. 
About ten basis states are sufficient to solve the Salpeter 
equation with rather high accuracy. The usual choice 
of 3-dimensional harmonic oscillator functions is less fa- 
vored here since their asymptotic behavior ~ e~ y ' for 
y — > 00 turns out to be not appropriate for deeply bound 
states like the pion. For other mesons, however, we have 
checked that both basis systems give equal results within 
numerical errors. 
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2. Expectation values 

I n o rder to compute the matrix elements present in 
cq.(B6) it is useful to rewrite (ip\il>) and (i/)\Hip) in terms 
of $^ and <f> using eq.(p2). We find 



(V#> = 



2lU\ LU 2 



cui m 2 + u) 2 mi 
tr (($++) t $ ++ - ($~)t*- 



(B8) 



for the norm. The matrix elements of H. can be split as 



(1>\H1>) = (1>\Til>) + <</W> 
(iplTip) = J (wi+«2)tr [$ f $] 

(V>|VV>) = - f f tr [$ t 7°(^$ , )7°] 



(B9) 
(BIO) 

(Bll) 



with J = J d 3 p/ (2ir) 3 . The basic formula for the cal- 
culation of the kinetic energy and confinement matrix 
elements is given by 

tr($ t $')= (B12) 
= tr[($++) t ($ ++ )' + ($— ) t ($— )' + 

+($+-)t($+-y + ($-+)t($-+)'] = 

= tr[($++) t ($ ++ )' + ($— ) t ($— )' + 

($++) t ($ ++ )' s's - ci4 ($++) t )'s - 

-c 2 ci ($— ) f s($++)' a' + c 2 c' 2 ($~)t ss '($— )' + 
-hcic; )t(*~ )' s's - ci4 ($~ ) f s'($++)'s - 
-020; ($++) t «(*—)' s' + c 2 c' 2 ($ ++ )t ss' ($++)'] 

where the prime indicates the dependence on p '. To com- 
pute the kinetic energy term eq.(BlC) set p' = p in this 
equation. 



For the interaction term eq.( Bll ) we investigate the 
following contributions: 



a. Scalar confinement 
From the scalar confining kernel eq.(||) one has 

&\V%i))=- [ f V c tr [$t/$' 7 o] (B13) 

with V c = Vc((p-p') 2 )- Since 

7»*7» =(_t:- (B14) 

the expression for tr [$t 7 7 ] is obtained by chang- 
ing the sign of the first two terms ($ ++ )t($++)' and 
($— )t($— )' in eq.(|l|). 



6. Vector confinement 
From the vector confining kernel eq.(^) one has 

(V>|V£V) = // V c tr [$t$'] (B15) 



and eq.(||T|) can be applied directly. Note that a covari- 
ant form to write the vector confinement kernel is 



[K(P,p,p') X (p')} 



(B16) 



with p± = p — (Pp)/P 2 , so that the relation 
P^[K(P,p,p')]=0 



(B17) 



holds as is required to rewrite the normalization condi- 
tion for <F 



c. 't Hooft interaction 

From the 't Hooft kernel eq.(^0|) we find (omitting fla- 
vor indices) 



(B18) 

tr (l 4 tr$' + 7 5 tr($'7 5 )) 7 ] = 

(tr$)* (tr$') - (tr$7 5 )* (tr$' 7 5 ) 



= -AG 

= -AG 
with 



Using the decomposition eq. (jB3j) and tr ^5 q 
one obtains 

tr <&++>-) = V2n ( +l (p)Yj M A^p)Sso 

tr ®jmj (f) = V2 1ZP J0 (p) Yj Mj (Sip) 5 S 

We further see that 

tr(cp Sq ap) = tr(apip Sq ) = V25siP q 



tr$ = tr($+~ + $~+) 
tr($7 5 ) = tr ($++ + $ ) 



(B19) 
V^Sso 

(B20) 
(B21) 

(B22) 



with the spherical components p q = \J An ji pYi q (fl p ) 
We find 



dO p tr($7 5 )= (B23) 
V2 VA^ (K% (p) + K { ~} (p) ) 6 S 

d% tr $ = (B24) 
-V2 VA^p (ci + pa) (n£} ( P ) - TZ { N J (p)) 5s 1 8 L 1 <5j 







3. Matrix elements 

The expressions obtained above can now be used to 
compute the required matrix elements. In order to calcu- 
late e.g. = (e- + ^ \bj) one has to replace in eq.([Bg|) 

($++) t by E\ and $++ by Ej setting ) f = and 
= 0. In the following we will use the notation 

(i\f(p)\j) = J J^Rn.LMRn jL Ap)f(p) (B25) 

for the radial integrals between the basis states. These 
integrals can be effectively computed numerically using 
e.g. Gauss quadrature routines. According to eq.(^) we 
write Hg = T?f + Vff with V t f = (V c )?/ + (V T )i( ■ 



Si2(L' , S , L, S, J) — 

= £ Si(L',5',L fcs 5 fc ,J) S 2 (L k ,S k ,L ) S ) J) (B31) 

L k S k 

The kinetic energy matrix elements are now given by 

rp~\ — h rp 

ij ~ V ~ 

= (i|( Wl + ua) (1 + p 2 (cj + c 2 2 j) \j) 5 LiL] 6 SiSj 

rp^ rp h 

ij ~ ij ~ 

= {i\(u 1 +uj 2 )2p 2 c 1 c 2 \j)S 12 (i,j) (B32) 
with S 12 (i,j) = S\ 2 (L l , Si, Lj,Sj, J). 

c. Scalar and vector confinement matrix elements 



a. Normalization matrix elements 
For the normalization matrix elements we find from 

2uj\uj 2 



eq.(gS 
N 



-N: 



uim 2 + uj 2 mi 







(B26) 
(B27) 



b. Kinetic energy matrix elements 

In order to evaluate the angular momentum structure 
of matrix elements the relation (—io~ 2 ) a (ia 2 ) — — a is 
quite useful. Define the qq spin matrix as xs q = fs <j 
and use s = op. Then one can write e.g. 

tr (<PSi Mi s fS, Mj S) 
= ~ tr (X<?» Mi Ol ® s 2) XSj Mj ) = 

=: - (SiMi \si <g> s 2 \ SjMj) (B28) 

with Si — 2Sip where <Si is the spin operator acting on 
the first quark and S 2 on the second quark. It is useful 
to define 

Si(L',S',L,S,J) := ([L'®S'} J \ Sl / P \[L®S} J ) (B29) 
S 2 (L', S', L, S, J) := ([L' ® S'] J \s 2 /p\ [L ® S] J ) 
S U (L', S', L, S, J) := ([L 1 ® S'] J \ {s x ® s 2 )/p 2 \ [L ® S] J ) 

Using the Wigner-Eckart theorem and the notation L = 
\/2L + 1 one has 

Si{L', S', L, S, J) = (-1) 5+5 ' S 2 (L', S', L, S, J) = 
= {-l) L ' +1 JL'LS'S V9V12 



V 1 L 






5" 






S 






1 





1/2 1/2 S' 
1/2 1/2 S 
1 1 



(B30) 



The confinement matrix elements can be computed by 
inserting two complete sets of basis functions like 

(i\A(p)V(r)f 2 (p')\j)= (B33) 
= £<*1/i(p)l<?> (9\V(r)\h) (h\f 2 (p')\j) 

where the matrix element of V can be evaluated in coor- 
dinate space. One finds for the scalar confinement matrix 
elements 

(vS)t, + = (yS)Tf = (i\v c \j) - £ ( B34 ) 

g,h 

{ (i\p Cl \g)S 2 (i,g) (g\V c \h) (h\p Cl \j) S 2 (h,j) 
+(i\pc 2 \g)S 1 (i >9 ) (g\V c \h) (h\pc 2 \j) Si(h, j) } 



and 



(VS)tj- = (VS)^ + = -Yl (B35) 



y.li 

{ (i\p Cl \g)S 2 {i,g) (g\V c \h) (h\pc 2 \j) S^h, j) 
+(i\pc 2 \g)Si(i,g) (g\V c \h) (h\p Cl \j) S 2 (h, j) } 

with 

(g\V c \h) = fr 2 drR? HLt (r)RF jLj (r) V£(r) 6 L(Lj 6 SiSj 
where R F denotes the Fourier transformed basis func- 
tions and V^(r) = a c + b c r. 

For the vector confinement matrix elements one only 
has to change the sign of h in the two equations 
above. 



d. 't Hooft matrix elements 

The 't Hooft matrix elements with the regularizing po- 
tential can be computed analogously to the confinement 
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matrix elements. In the unregularized case A — > cor- 
responding to Vf e g(r) — > S(r) only L = basis states 
contribute to (g|Vreg|^)- Consistently only L = states 
will be taken into account for (g|Vreg|^) also for A > 
so that the angular selection rules are not changed by the 
regularization. The result for L = S = reads 

(V^ + )f=° = (F T -)r° = (V T + -)ff° = (V^ + )f=° 
= 8G (i|V reg |:0 S Si o S Sj o 6 Lt o ^ o <5jo (B36) 

and for L = S 1 = 1 

= -8G £ (ci + c 2 )\g) ( 5 |Vreg|/i)(/i|p (ci + c 2 )|i) 

5siiSSji5L t i8Lji8jo (B37) 

This result shows that the 't Hooft interactions affects 
only mesons with J = and L = S = (i.e. pseu- 
doscalar mesons with J Vp = 0~) or L = 5 = 1 (i.e. 
scalar mesons with J Tp = + ). In the nonrelativistic 
limit the contributions to the scalar mesons vanish. 
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1 2 1 2 

FIG. 2. Instanton induced interaction vertices for qq — > qq. 

FIG. 3. Mass spectra of the pseudoscalar mesons. The 
columns for each meson correspond (from the left) to model 
VI, model V2, experiment (pBJ and model S. The shaded ar- 
eas (3rd column) indicate the experimental full width of the 
meson. 

FIG. 4. Mass spectra of the vector mesons (see also caption 
to Fig. g). 
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FIG. 5. Regge trajectory for the isovector mesons with 
5 = 1. The solid line shows the experimental masses for 
p, a,2, P3, ei4 (5^] where the errorbar gives the experimental 
error for the resonance position. The short dashed line corre- 
sponds to the calculated masses of model VI, the dotted line 
to model V2 and the long dashed line to the scalar confine- 
ment (model S). 

FIG. 6. Effective mass dependence for the non strange 
(solid curve) and strange quark (dashed curve) on the in- 
stanton size cutoff p c . 



TABLE III. Parameters of the different models (see 
Sec. |lVA[ ) 



Parameter 


V 1 


Vz 




B 


m„[MeV] 


170 


340 


340 


m s [MeV] 


390 


568 


487 


a c [MeV] 


-552 


-1340 


-998 


b c [MeV/fm] 


570 


1400 


1000 


g [MeVfm 3 ] 


51.67 


34.65 


44.79 


g' [MeV fin 3 ] 


46.92 


30.84 


41.01 


A [fm] 


0.42 


0.42 


0.42 



FIG. 7. Effective coupling constants for the instanton in- 
duced non-strange quark interaction g (solid curve) and the 
strange non-strange quark interaction g' (dashed curve) as a 
function of the instanton size cutoff p c , see eq.< 



FIG. 8. Radial Pion amplitudes piz[^\p) (upper compo- 

dashed 



and p1Zq \p) (lower component, 



nent, solid curve 
curve) with the parameters of model VI 



TABLE IV. Spin orbit splitting for the first positive parity 
mesons and the effect of the instanton induced interaction on 
the ++ mesons, where a * denotes the existence of an addi- 



tional state with imaginary mass and zero norm, see Sec.tVB 
(all Masses in MeV). 



Meson 


J PC I 


VI 


V2 


S 


a (980) 


0++ 1 


960 


1130 


1260 


/o(975) 


0++ 


950* 


1270* 


950 


A (1285) 


1++ 


930 


1060 


1150 


/ 2 (1270) 


2++ 


1100 


1280 


1010 



FIG 



9. Radial Rho amplitudes pTZ\^ (j>) (upper 



component, solid curve), pTZ^ijp) (lower s-wave component, 
p7?-2i (p) (upper d-wave component, 



long dashed curve) , 
dashed dotted curve 
short dashed curve) with the parameters of model VI 



pTZ^i ] (p) (l ower d-wave component, 



TABLE I. Flavor dependence of the instanton induced in- 
teraction G/ 2 f 3 j 1 f i (see eq.(p2|)) 



hh -> 

/a/4 I 


ud du us 


ds 


sd 


s u 


uu dd ss 


ud 
du 


-g 
-g 










us 
ds 
sd 
su 


-g' 







-g' 







-g' 








-g' 




uu 
dd 
ss 










g g' 
g g- 
g' g' 


TABLE II. Flavor matrix elements of Gf 2 f 3 j 1 f i for pseu- 
doscalar mesons 




n Q 




T]n 






7r u 

'In 

Vs 


-g 








g 

V2g 


1 




^2g' 




TABLE V. Comparison of experimental and calculated 
meson decay observables for the Salpeter models VI, V2, S 
and nonrelativistic results NR 



Mesonic decay experimental [25 VI V2 S NR 



U [MeV] 
Sk [MeV] 

r(7T° -» 77 ) [ 

T(r) -» 77) 
T(r)' -> 77) 
r(p -> e + e~) 
r(uj -» e+e~] 



[eV] 
eV] 
[eV] 
[keV] 
[keV] 
[keV] 



131. 7± 0.2 
160.6 ± 1.4 
7.8 ± 0.5 
460 ± 5 
4510 ± 260 
6.8 ± 0.3 
0.60 ± 0.02 
1.37 ± 0.05 



130 260 

180 300 

7.6 4.0 

440 220 

2900 2030 

6.8 28 

0.73 3.1 

1.24 4.5 



200 
210 
4.4 
220 
1390 
8.1 
0.87 
1.50 



1440 
730 
30000 
18500 
750 
8.95 
0.96 
2.06 



TABLE VI. -q, 7] mixing parameters from BS norm (see 
cc l- (JH)) compared to data calculated from experimental J/'I' 
decays fcol 



Meson 


mixing coefficient 


J/^ decay 


VI 


V2 


S 


r;(547) 


1^1 


0.63±0.06 


0.71 


0.71 


0.70 




inl 


0.83±0.13 


0.70 


0.70 


0.72 


V(958) 


|AVI 


0.36±0.05 


0.85 


0.78 


0.83 




IV 1 


0.72±0.12 


0.52 


0.63 


0.55 
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